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Abstract. The paper focuses on a model describing the spreading of an insoluble surfactant on a thin viscous film 
with capillary effects taken into account. The governing equation for the film height is degenerate parabolic of fourth 
order and coupled to a second order parabolic equation for the surfactant concentration. It is shown that nonncgativc 
weak solutions exist under natural assumptions on the surface tension coefficient. 



1. Introduction 

The modeling of the spreading of an insoluble surfactant on a thin viscous film leads to a coupled system of 
degenerate parabolic equations describing the space and time evolution of the height ft. > of the film and the 
surface concentration T > of surfactant. Assuming the film thickness to be small enough so that lubrication theory 
is applicable, taking into account capillary effects but neglecting gravitational and intermolecular (van der Waals) 
forces, the following system is obtained [H El IE] 

9th + d x Q h 3 dlh + i h 2 d x a{Y)\ = 0, (t, x) £ (0, oo) x (0, 1) , (1) 

d t T + d x Q h 2 r dlh + h r d x <r(T) ) = Dd 2 r, (t, x ) e (o,oo) x (0,1), (2) 

with homogeneous Neumann boundary conditions 

d x h(t,x) = d x h(t,x) = d x V(t,x) = 0, (t, x) G (0, oo) x {0, 1} , (3) 

and initial conditions 

(ft, r)(o) = (/*>, r ), xe(o,i). (4) 

Here, <j(T) denotes the surface tension which depends on the local concentration of surfactant, and D > stands for 
the surface diffusivity of the surfactant. Since a is defined up to a constant, we may assume without loss of generality 
that cr(l) = as a normalization condition. As the presence of surfactant reduces surface tension, a is a non-increasing 
function of T; for instance, 



a (s) := 09 + 1) 



1 - s 



+ l^ 1/3 



-13, s>0, (5) 



P 

with /3 € (0,oo) and its limit as /3 — > oo (which is often assumed in applications) 

a-oo(s) := l-s, s>0. (6) 

The system (jTJ)- (J3J) is a fully coupled nonlinear system of parabolic equations featuring a degeneracy where h vanishes, 
a fact which cannot be excluded a priori. Thus, classical solutions are unlikely to exist for all times in general and only 
local existence of smooth solutions to (HJ-© in the absence of capillarity have been shown in [TUJ[TT]. The alternative 
is to study the Cauchy problem in a framework of weak solutions (see Section [2] for a precise definition) and this 



This work was partially supported by the french-german PROCOPE project 20190SE. 



2 



JOACHIM ESCHER, MATTHIEU HILLAIRET, PHILIPPE LAURENCOT, AND CHRISTOPH WALKER 



approach has been successfully employed to establish the existence of weak solutions to systems similar to ([T])-©, for 
instance when the capillarity effect are neglected but the gravitational ones are accounted for [S] or when Y is replaced 
by A(r) = max{0, 1 — T} + 1 in © [TJ|2], a finite element numerical scheme being also developed in these two papers. 

To our knowledge, existence of weak solutions has only been tackled in [7J where this is proved under further 
technical assumptions on the surface tension a. Namely, given a surface tension a £ C ; 2 o *(R), defining the free energy 
9a by 

g a {l) = g' a (l) = 0, g»(s) = for s G R, (7) 

the authors assume the following: 
(A4) The function g a lies in C?£(R). 

(A5) There exists c g > such that g'l(s) > c g for all s £ R. 

(A6) There exist C g and some r £ (0, 2) for which g£(s) < C g (\s\ r + 1) for all s G R. 
The need in [7J to define a and g a in R instead of the physically relevant range [0,oo) for surfactant concentration 
stems from the fact that the weak solution (h,T) to ([I])-© constructed in [7J might not satisfy Y > 0. The extension 
of a and g a to negative values then induces several limiting conditions. Namely, the convexity (A5) of g a in R requires 
not only that cr'(s) < for s > as expected but also cr'(s) > for s < 0, which implies cr'(O) = and thus excludes 
surface tensions a like ap in ([5]) and CToo in (|6]). In addition, assumption (A5) yields cr'(s) < — c g s for s G R so that a 
necessarily has a quadratic decay at infinity. This again excludes the previous examples ap in ([S]) and (Too in ([B]). 

The aim of the present paper is to construct a weak solution to (H])-© under weaker assumptions on the surface 
tension a (satisfied in particular by a^) and such that both h and Y arc nonncgative throughout time evolution. More 
precisely, we assume that the surface tension a satisfies: 

(HI) a G C 1 ((0,oo)) nC([0,oo)) with cr(l) = 0. 

(H2) There exist cr , <ti G (0, oo) and 9 £ [0, 1) for which: 

-(T < o-'(s) < -— ^ fors>l, -o- <a'(s)<0 for se (0,1). (8) 
1 ^ s 

The assumptions (H1)-(H2) include physically relevant surface tensions er, which may slowly decrease at infinity. In 
particular, Ooo is included (but not 073 for /3 G (0,oo)). 

In the next section, we introduce the definition of weak solutions and give precise statements for our existence result. 
To prove this result, we first construct nonnegative solutions in the framework of [7J under assumptions (A4)-(A6). 
This improves the results of [7j in that the surfactant concentration Y stays nonncgative through time evolution. This 
is the content of Section [3] Then we extend the construction to a surface tension a merely satisfying (H1)-(H2) by 
approximating a with surface tensions satisfying (A4)-(A6) and studying compactness properties of their associated 
weak solutions. The construction of er^ and the compactness argument are presented in Section [4] 



2. Weak solutions and main results 

To introduce the definition of weak solutions, we first derive energy estimates satisfied by smooth nonnegative 
solutions to (HJ) . So, let us consider a non-increasing and smooth surface tension a and a smooth solution {h, Y) to 
(JT])-(j3l) in (0,T) x (0, 1) for some T > 0, both functions being uniformly bounded from below by a positive constant. 
First, we note that by (JT|)- ([3]) there holds 

- [ \\ d* 

Setting 

fc,(l)= 5 ;(l) = 0, g'l(s)=-^M. for S G(0,oo), (10) 



d 
dt 



r dx 



(9) 
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it follows from (HJ)-© that 
d 
dt 



\d x h\ 2 



+ g a {T)\ dx 



d 2 x hd x 
h 2 



h 3 h 2 
^d x h+^d x a(T) 



dx 



r d 3 h + h r d x a(T) - D d x T 



dx 



1 r /? 3 a'(V) 

- \d 3 x h\ 2 + h 2 d x a(T) dlh + h\d x a(T)\ 2 - D-^\d x T\ 2 



n 



r 



dx. 



We abbreviate product terms by introducing J 2 and J 2 , where 

h 3 ' 2 



This yields 



d 
dt 



J f = J f [h,T] : = 
J s = J s [h,T] := 

\d x h\ 2 



3 

h 3 ' 2 



h 1 / 2 

3 x h+— d x <j(T), 



dlh + h 1 / 2 d x a(T). 



(11) 
(12) 



+ g CT (r) dx 



l\J f [KT]\ 2 + ^\J s [h,Tf 



h 3 \dlW 2 



24 



1 



h\d x a(T)f 



dx 



(13) 



D 



V(T) 



\d x r\ 2 dx. 



Consequently, we infer that, regardless the qualitative properties of a, (h,T) should satisfy 

fteM0,T; ^(0,1)), re 2^(0, T;L X (0,1)) 

together with 

h 3 ' 2 dlh e L 2 ((0, T) x (0, 1)), h l / 2 d x a{T) e L 2 ((0, T) x (0, 1)). 



(14) 
(15) 

Since D > 0, the energy estimate (|T3|) provides an additional estimate which depends strongly on the properties of a, 
namely -J— a'(T)/Td x T 6 L2((0,T) x (0, 1)). We will actually prove that, under assumption © on cr, this additional 
estimate guarantees that the solutions we construct satisfy the further regularity 

rei 2 ((0,T)x(0,l)) and a(T) € L 4 / 3 (0,T; W* /3 (0, 1)), (16) 

see Lemma [9] Let us point out here that © implies that er' docs not decay too fast towards -oo at infinity. 
Hence, assuming ([14"]) and (Ti~5j) to hold true, we realize that Jf[h, T], J s [/i, T] £ L 2 ((0 7 T) x (0, 1)). Since an alternative 
formulation of CD)-© reads, 

d t h + d x (h 3/2 J f [h,T]j = in (0, oo) x (0, 1) , (17) 

a,ri4(/i I/2 rj s [)i,r]) = d d 2 x v in (o, oo) x (o, 1) , (is) 

we infer from i[H j) -([ll) ]) and the embedding of iJ^O, 1) in 1^(0, 1) that h 3 ' 2 J f [h, T] and h 1 / 2 T J s [h, T] both belong- 
to Li((Q,T) x (0, 1)) and we can give a meaning to ©-© at least in the following weak sense: 

Definition 1. Let T > and a be a surface tension such that either 
• a e C 1 (0,oo)nC([0,oo)), cr(l) = 0, and © holds true, 

or 
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• a G C 1 ' 1 (R) is such that er(l) = and g a , defined in ([7|) 7 satisfies (A4)-(A6). 

Then, given an initial condition (ho, To) G ^ 1 (0, 1) x £2(0, 1) with ho > and To > we say that (h,T) is a weak 
solution in (0,T) to |T])-(j4|) with surface tension a and initial condition (ho, To), if 

• h > and T > satisfy 

h G Loo(0, T; if^O, 1)) n C([0, T] x [0, 1]), T G L^O, T; L^O, 1)) n L 2 ((0, T) x (0, 1)), 

d x h <E L 2 (V h (S)) for all 5 > 0, <r(T) G L x (0, T; 1^(0, 1)), (19) 

^^fe g L 2 (n), /^^(r) G L 2 ((0,T) x (0,1)), 

where Vh(5) '■= {{t, x) G (0, T) x (0, 1) : x) > 6} for S > andVh ■= {(t,x) G (0, T) x (0, 1) : h(t,x)>0}. 

• /or any ( € C°°([0,T] x [0,1]) swc/i t/iai C(T,a;) = for all x G [0,1] and d x ((t,x) = /or afl G 
[0,T] x {0,1}, there holds: 



T r l 



JO 



hd t c 



i h 3 1(0,00)^)^ + 5 /i 2 a x a(r) 



9 X C^ dxds = — y /io(a;)C(0, x)da; . 



(20) 



and 



lo lo \ 



- /i 2 l 



(0,oo 



(/i)r a^r) 



(x)C(0,x)dx. (21) 



With these conventions, our main result reads: 



Theorem 2. Lef £/ie surface tension a G C (0,oo) n C([0, oo)) satisfy o(T) = and (|8j). Then, given an initial 
condition {h ,T ) G ^(0,1) x L 2 (0, 1) with h > 0, T > and any T > 0, there exists at least one weak solution 
(h,T) in (0,T) to (HJ) — d3J with surface tension a and initial condition (ho, To) in the sense of Definition^ 

As mentioned in the Introduction, we split the proof of Theorem [2] into two parts. First, we focus on the 
nonnegativity issue of solutions to (H}-©. In this respect, we go back to the framework considered in [7] and we 
prove: 

Theorem 3. Let the surface tension a G C 2 (K) be such that c(l) = and the free energy g a defined by (JT)) satisfies 
(A4)-(A6). Then, given an initial condition (ho, To) G 7? 1 ((0, 1)) x £ 2 ((0, 1)) with ho > 0, To > 0, and any T > 0, 
there exists at least one weak solution (h,T) in (0,T) to (JT)) — with surface tension a and initial condition (ho, To) 
in the sense of Definition^ 

Moreover, the solution satisfies the further regularity 



T G Loo(0,T; L 2 (0, 1)) H L 2 (0,T; H\0, 1)) 



and the energy estimate 



where 



sup 

*e[o,r] LJo 



\d x h(t,x)f 



+ 9a(T(t,x)) 



dx\+T>[h,T] < j 



\d x ho(x){< 



+ 9a(T (x)) dx 



(22) 



(23) 



V[h,T] := ££ ( ^^^f ^ X)) \ \8Mr, X )\ 2 + h ^\dMnr, X ))\^ «fedr 



-D 



o Jo 



i „/ 



a'jT^x)) 
T(t,x) 



\8 x T(t,x)\ 2 dxdr. 
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With the regularity (|22|) . any weak solution constructed in [7] is a weak solution in our sense (see the proof of 
Theorem [3] for further details). The major novelty in this result is that we obtain nonnegativity of the surfactant 
concentration T. For the proof of Theorem [2j we consider a surface tension a £ C([0, oo)) n C 1 (0, oo) and introduce 
a family of approximate surface tensions (cr k ) k£ ^ satisfying the assumptions of Theorem [3] We achieve our result 
by studying the compactness properties of the family of associated weak solutions. A fundamental argument will be 
that, owing to assumption ((5]) and equation (|13p . the dissipation of energy is measured by 

ldxTl \dxdT. (24) 







r(i + 

When 9 £ [0, 1), this quantity enables us to control y/T in some Holder space (sec Lemma [9|). This, in turn, yields 
compactness on the concentration for any bounded family of solutions in ^((O, T) x (0, 1)). 

3. Existence of nonnegative solutions for decaying surface tensions 

In this section, we assume a £ C 2 (R) is such that tr(l) = and the free energy g a , as defined in ([7|), satisfies 
(A4)-(A6) and we construct nonnegative weak solutions to ([I])-©. In Sect. 3. 4], the authors remark that, for 
proving nonnegativity of the surfactant concentration of weak solutions to ([I])-((4]), a difficulty arises when multiplying 
equation @ by T_ = — min{0, T}. Indeed, under assumptions (A4)-(A6), the very low regularity of T implies only 
that d t T £ L 3/2 (0,T; (W£(0, 1))*) and T_ £ L 2 (0, T; H^O, 1)). This regularity docs not allow to define the duality 
bracket (dtT, To construct weak solutions with nonnegative surfactant concentrations, we go back to the strategy 
applied in [7] : construction of solutions to a regularized problem via a Galcrkin method, followed by a compactness 
argument when the regularization parameter goes to 0. We introduce a supplementary truncation operator in the 
regularized problem in order to guarantee that the solutions to the regularized problems have nonnegative surfactant 
concentrations. 

Throughout this section, we fix a nonnegative initial condition (ho, To) £ H (0, 1) x £2(0, 1). We also introduce a 
Lipschitz continuous truncation function T such that 

( s if s £ (0,1), 

T(s)=i 2-8 ifse[l,2], T(-s) = -T(s) if s < 0, (25) 

[0 if s > 2, 

and put 7fc := kT(-/k) for k > 1. Then, we set 

:= / T k (a'(r))dr for s £ R. (26) 



Ji 

We emphasize that this construction ensures that £ C la (R) has bounded first and second derivatives. Associated 
to this truncation of <r, we introduce a truncation of the identity 

Tk(s) := s ^4 for s e M - ( 27 ) 
We note that the construction above is well-defined because 

> oi(s) > o-'(s) for all s £ R. (28) 
With these conventions, our regularized problem reads 

d t h + d x ([a 3 (h) + l/k] d 3 c h + a 2 (h) d x a k (T)) = 0, (t, x) £ (0, 00) X (0, 1) , (29) 

d t T + d x (a 2 (h) T k (T) d 3 x h + ai (h)T d x a k (T)) = D d x T , (f, x) £ (0, 00) X (0, 1) , (30) 

subject to ©-(HI), where k is a positive integer. The notation a,(/i) stands for (max{0, h}) z /i for i = 1, 2, 3. This is 
the same convention as in [7] so that conditions (A1)-(A8) therein are satisfied. 
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3.1. Existence for (|29|l - (|30|) . To begin with, we fix k > 1 and prove: 

Lemma 4. Consider an initial condition (ho, To) G ^((O, 1)) x L2((0, 1)) with ho > and To > 0. For any k > 1 
and T > 0, there exists at least a couple of functions (h,T) having the regularity 

h e Loo(0, T; H\0, 1)) n L 2 (0, T; ff 3 (0, 1)) , T g ^(0, T; L 2 (0, 1)) n L 2 (0, T; H^O, 1)), (31) 
d^eL^T;^ 1 ^,!))*) , ^^3/2(0,^(^(0,1))*), (32) 

and satisfying, 



(d t h, ()ds- / (a 2 (h)8 x a k (T) + [a 3 (h) + l/k] d 3 x h) d x ( dxds = 0, 
Jo Jo 

for all C € L 2 (0,T;H 1 (0, 1)), together with 

[ (d t T,Qds- [ [ { ai {h)T d x a k (T) + a 2 (h) r k (T) d 3 x h- Dd x T) d x ( dxds = 0, 
Jo Jo Jo 

for all C G L 3 (0,T;W£(0,l)) and 

(h(0,-),T(0,-)) = (h o ,T ), 

the latter being meaningful as h € C([0,T]; (H l (Q, 1))*) and T E C([0, T]; (^^(0, 1))*) by flST]) and ([521) . 
Moreover, there holds the energy inequality 



sup 

te[o,T] 



l r 



\d x h{t,x)f 



g a (T(t,x)) 



dx } +V k [h,T) < 



\d x ho(x)f 



ga{T {x)) 



dx, 



(33) 

(34) 
(35) 

(36) 



where 



f> k [h,T} := 





"1 03(h)' 




k 7 



\d 3 x h\ 



D^\d x T\> 



+ ^\d x a k (T)f 



dxds. 



Remark 5. Note that, in (|36[) . <j k only appears in the last term of T> k [h,T\. 



Proof. We follow here the Galerkm method from [71 Section 3]. The system ([21?)) - ([3T))) is actually almost identical 
to the regularized system used in [7J Section 3] except that the truncation function r k is replaced by the identity 
there. Since r k is a bounded and Lipschitz continuous function, the analysis performed in [Jj Section 3] carries over to 
([29]) - (j30|) with only slight changes, the main one arising in the derivation of the energy inequality. We will thus only 
give a sketch of the proof and refer to [7J Section 3] for details. 

The first step is an alternative formulation of ([2^[) - ([5n[) in therms of h and the new unknown function v := g' a (T), the 
latter being well-defined thanks to the convexity (A5) of g a . Denoting the inverse function of g' a by W, we have 

d t h + d x ([a 3 (h) + l/k] d x h-a 2 (h)T k (W(v))d x v) = 0, (37) 

d t W(v) + d x (02(h) r k (W(v)) 8 3 x h ~ oi(fc) W(v)r k (W(v))d x v) = D d 2 x W(v) , (38) 

in (0,oo) x (0, 1). As already mentioned, ([3"Tj) - ([3"5)l is the same as the system studied in [7J Section 3] except for the 
terms involving the bounded and Lipschitz continuous function r k . Not surprisingly, considering the same Galerkin 
approximation to ([37 [) - (j38)) as in [3 Section 3], one can prove the local existence of solutions to the Galerkin approx- 
imations exactly in the same way as in [71 Section 3.1]. To obtain the global existence, we argue as in [7J Section 3.2] 
by deriving an energy estimate for the Galerkin approximations. Since there is a slight modification necessary, let us 
sketch the proof for ([3"7[) - ([3"5)) . the argument being the same at the level of the Galerkin approximations. We multiply 
(|3"7[) by —d x h, (|38[) by v — g' a (T), integrate over (0, 1), and add the resulting identities to obtain 



dt 



\d x h\'< 



9a{T) 



Ax 



^ +a 3 (h) 



\d 3 x h\ 2 - D ^Pl^lf + ai (h)(a'(T)a' k (T))\d x Tf }dx = I, 
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where (see ([7]) and ([2? 



03(h) 8 X T d 3 x h(a' k (T) + a'(T) 



dx 



2a 2 (h) d x a k (T) d 3 hdx. 



Since a' k < 0, it follows from ([25} that (o~' k ) 2 < c'oj. while Young's inequality ensures that 

\2a 2 (h) d x a k (T) d 3 h\ < \d x a k (T)\ 2 + ^1 \% h \* , 



so that we finally obtain 



dx 



7 + a 3 (h) 



\d 3 x h\ 2 - D^p-\d x T\ 2 + ai (h)\d x a k (T)\ 2 } dx 



r 



< 



da; . 



This yields 



\d x h\ s 



dx + 



a 3 (h) 



3 t |2 



D^l\d x r\ 2 + ^ |cwr)| 2 } dx < o, 



dt J [ 2 
whence (|36j) after time integration. 

The convergence of the Galcrkin approximations to a solution to ([29 |) -([30 l) satisfying the properties listed in Lemma |4] 
is then carried out as in [7J Section 3.3] to which we refer. □ 

At this point, we show that the idea to introduce truncation functions T k and o~ k yields the nonnegativity of T. This 
relies on a gain of regularity for <9 4 r. 

Lemma 6. Consider an initial condition (ho, To) G if 1 ((0, 1)) x L 2 ((0, 1)) with ho > and To > 0. Given k > 1 and 
T > 0, any solution (h,T) to ([3l])-([36]) and ©-© m fee sense 0/ Lemma @] satisfies d t T G L 2 (0,T; (if^O, 1))*) and 
T > a.e. m (0,T) x (0, 1). 

Proof. Owing to (f3~Tj) . the embedding of i? 1 (0, 1) in Lqo(0, 1), and the compactness of the supports of a' k and r k (which 
follows from (A5) and the properties of T), there holds 

a,(h) T O x a k (T) + a 2 (h) r k (T) d 3 h = ai (h)Ta' k (T)d x T + a 2 (h) r k (T) d 3 x h G L 2 ((0,T) x (0,1)), 

and Dd x T G L 2 ((0,T) x (0,1)). As a consequence ((HI) also holds true for all C G i 2 (0, T; iJ x (0, 1)) and <9 t r G 
L 2 (0,T;(i? 1 (0,l))*). Then, if /? G C 2 (M) is such that fi' is Lipschitz continuous, we have (i'(T) G L 2 (0, T; if^O, 1)) 
and 

f /3(r)da; = (a t r,/3'(r)). 



dt Jo 

Assuming furthermore that (3 is convex, i.e. j3" > 0, there holds, for any t G (0,T), 



/3(r(i))dx< / p(T )dx + 







(01(h) r d x a k (T) + a 2 (h) r k (T) d 3 x h) p"(T)d x T\ dxds. 



To finish off the proof, we apply this inequality to a family of functions approximating the negative part of T. Namely, 
we fix a nonnegative \ G Cq 



'S°( 



such that x ^ has support in ( — 1,0) and define /?i by 
XT X(") d " 



A(0)=0, 



for s 6 
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We then set j3 e (s) :— £/3i(s/e) for s € K and e > 0. Taking j3 — f3 e for e > in the above inequality, there holds, for 
each t g (0,T), 



&(r(i))da: < 



T i-l 



^0 



(oi(/i) r d x a k {T) + a 2 (h) Tfc(r) ^/i) /^(iWl dxds : 



since /3 e (r ) = due to T > 0. Observing that \T k (s)/3' s '{s)\ < \sP' e '(s)\ < C(x) for s € R, we have 



/ / | (oi(ft) r d x a k {v) + a 2 (h) rfc(r) a^) ^'(D^rl dxds 

Jo Jo 



<C(x) 



{|r|<e} 



|ai(/i) fllcfffc (r)9 a r| + \a 2 (h) dlh d x T\] dxds 



<C( X ) [||oi(/i)a / fc (r)|| Loo( o > r ;ioo (0,l)) + ||a 2 (/ l )|U ao( o 1 T;L ao (0 > l))] 



{|r|<e} 



|^r| 2 + \dlh d x v\ 



dxds. 



As d x T and d x h both belong to L 2 ((0,T) x (0, 1)) and d x T = a.e. in {V = 0} by [§J Lemma A.4], we obtain in the 
limit e — > 



/ max{-r(i,x),0}dx < 0, for ie(0,T). 
Jo 



This completes the proof. 



□ 



3.2. Proof of Theorem[3l Let a be as in the statement of Theorem [3] and consider an initial condition (ho, Tq) g 
i? 1 ((0,l)) x L 2 ((0,1)) with /i > 0, L > and T > 0. First, applying Lemma [4] and Lemma [6l we obtain a 
sequence (/i fe ,r fc ) fc >i of solutions to p ^l -pO ]) . ©-© for which d t T k g L 2 (0, T; (£^(0, 1))*) and T fc > a.e. in 
(0,T) x (0, 1). In particular, for each k > 1, the time regularity of h k and F^, together with the initial conditions ff35|) 
(hk(0, -),T k (0, •)) = (/lojTo), yield the integration by parts formula: 



(dth k ,0 dt 



(d t T k ,0 dt 



h (x)({0,x)dx 



r (x)C(o,x)dx 



T r l 



h k (s,x)d t ((s,x)dxdt , 



o Jo 

T r l 



JO 



T k (s,x)dtC(s,x)dxdt 



for any test function C g C°°([0,T] x [0, 1]) such that C(T,x) = for all x g [0, 1] and d x ((t,x) = for all G 
[0,T] x {0, 1}. Hence, taking such a test function £ in (|3"3")) - (fM)) we obtain : 



t r i 



o Jo 



a*3(h k ) + -j-) dlh k + a 2 (h k ) d x a k (T k ) 



d x ( ) dxdt 



/io(x)£(0, x)dx . 



(39) 



T /■! 



Jo 



(r fe OtC+ [a 2 (h k ) r k (T k ) dlh k + oi(hfc) T fc 9 xCTfe (r fe )] 9 x C + £T fe S|C) dxdi 



r o (x)C(0,x)dx. (40) 



So, the proof reduces to find a weak cluster point (h, T) of the sequence ((h k ,T k )) k >i that has the regularity (fT9|) and 
for which we can pass to the limit in the two previous equations. 

First, we note that the conservation laws © are also satisfied by (h k ,T k ). Consequently, due to (|3l))) and the 
Poincarc inequality, we have uniform bounds for 

• (fe fc) fc >i in Loo(0,T; # 1(0,1)) and {ga{T k ))k>i in £ qq(0, T; Xi(0, 1)), 

• (v^N^^jtM, (v^t^fc)SxO*(Tfc))fc>i, and (V"^^)/^^*)^! in £ 2 ((0,T) x (0,1)). 
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Owing to the bound (A5) from below on a', this yields a uniform bound on (Ti-)k>i hi Loo(0, T; £2(0, 1)) and 
£2(0, T; i? 1 (0, 1)), and the sequence of fluxes, given by 



t k 



a>2{hh) T k (T k ) 3 1/2 
77 u/2 — p — dZ,h k + ai{hk) 1 d x a k (T k ), 
ai(h k y/ z T k 



a 3 (h k 



1/2 



— ) d 6 x h k + a 2 (h k ) [ 3a 3 {h k ) + - 



-1/2 



d x cr k (T k ), 



are also bounded in L 2 ((0,T) x (0, 1)) by fl3gj) . 

Repeating the arguments in [3J Section 2] and Section 3.4], we may extract a subsequence (not relabeled) and 
find functions ft and T such that the following convergences hold: 

• h k ->hmC{[0,T ] x [0,1]) and T k ->• T in L 2 (0,T; L p (0, 1)) for all p£ [l,oo), 

• v/3ffl3(ftfc) + 3/kd*h k -± H in L 2 ((0,T) x (0, 1)) with # = y/3a 3 {h)d*h a.e. in {ft ^ 0}, 

• d x T k ^ d x T in L 2 {{0,T) x (0,1)). 

Arguing as in the proof of [7J Equation (3.28)], the previous convergences imply that 

(3a 3 (h k ) + 3/k)d x h k ^h^ 2 l (0tOs) (h)d x h in L a ((0,T) x (0,1)) . 

Next, interpolating the bounds on (T k ) k >i with the help of [5j Proposition 1.3.2], we deduce that (T k ) k >i is bounded 
in L 6 ((0,T) x (0, 1)) and that the convergence of (T k )k>i to T takes actually place in L p ((0, T) x (0, 1)) for allp e [2, 6). 
Now, since 

0<r fc (s)<s for all s > and r k {s) = s for < s < s k := \(k/C g ) r/{r+ ^ 
(see assumption (A6) and (|2"7[) for the definitions of r and Tfe, respectively), we have, for p > 1, 



1 



T 

Jo 



r fc (r fc ) 



dxdt = 



< 



{r fc > Sfc } 



r fc (r fc ) 



2'' 



{r fc > Sfc } 



- 1 



It dxdt < 



dxdt < 2 p 
C(p,T) 



dxdt 



{r k > Sk } 



Since s k — > 00 as k — > 00, we conclude that T k (T k )/T k — > 1 in L p ((0,T) x (0,1)) for any p > 1. Similarly, since 
cr^s) = cr'(s) for s e [0, s k ] and tr' g C j (K), it follows from (A6) and $Etfy that, given p a g [l,6/(r + 1)), i? > 1, and 
> 1 such that s k > R, we have 



< 



f fW k {T k )-a'{Y)r dxdt 
Jo Jo 

\a'(r k )-a'(r)\ P0 dxdt + 



'{max{r fc ,r}<ii} 



'{r fc >_R}u{r>_R} 

|r fc - rf° dxdt 



\a' k (T k )-a'(T)\ P0 dxdt 



C(p ,C gi r) 



{max{r fc ,r}<s fc } 



{r fc >_R}u{r>_R} 

T 



(r+l)p 



p(r+l)poj 

< m M f f\T k -vrdxdt^ c{p ^ r) 

Jo Jo 

< \\o"\\ Leel0 ,R) f T fn-rr dxdt + 



dxdt 



R e-(r+i) P0 J {Tk>R]u{r>R} 
C7( Po ,C7 3 ,r,r) 

^6-(r+l)p ' 



r b ) dxdt 
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Letting first k -> oo and then R -> oo yield that a' k (T k ) -> er'(r) in L Po ((0,T) x (0, 1)) for any p G [l,6/(r + 1)). As 
r < 2 we note that we may choose po > 2 in the previous convergence which, combined with the weak convergence of 
( d x r fc ) fc > i in L 2 ((0,T) x (0, 1) ) implies that d x ak(T k ) ^(r) in L go ((0,T) x (0,1)) for some g > 1- Consequently, 
y/a>i(hk) d x a k (T k ) y/ ' ai(h)d x cr(T) and 

-1/2 



a 2 (/ifc) 3a 3 (/ifc) 



9^ fc (r fe ; 



\/ai(h) 



d x <r(T) 



L 2 ((0,T)x (0,1)). 



Thus, we conclude that J s fc H/2 + y/ai(h)d x a{T) and J k f H + (y/a 1 {h)/2)d x a{T) in L 2 ((0,T) x (0, 1)). Combining 
these convergences with the convergence of (h k )k>i to h in C([0, T] x [0, 1]) and that of (Tk)k>i to T in L 2 ((0, T) x (0, 1)) 
allows us to pass to the limit in (|39 |) -(|40 |) . 

That h is nonnegative can be obtained as in Section 3.4] while the nonnegativity of T is preserved by the weak 
limit. Concerning the energy estimate (f23|) . we recall (A5) and prove as above that 



y/-a' k (T k )/T k -+ ^-o'(T)/T in L 2 ((0, T) x (0, 1)). 



Consequently, (\/—a' k (r k )/T k d x T k ) k >i converges weakly in Li((0,T) x (0, 1)) to y/— cr'(T)/Td x T, and we can pass to 
the weak limit in the energy estimate. This completes the proof of Theorem [3] 

4. Existence of nonnegative solutions for slowly decaying surface tension 

From Theorem [3] we obtain existence of weak solutions for a class of surface tension a decreasing at least quadrat- 
ically to — oo. We now extend with Theorem [2] the existence result to a class containing surface tensions which 
decrease slowly to — oo at infinity (but not too slowly, see (H2)) and are thus closer to applications. To this end, we 
fix a surface tension a satisfying (H1)-(H2) and an initial condition (ho, To) € i? 1 (0, 1) x i 2 (0, 1) satisfying ho > 
and r > 0. We split the proof of Theorem [5] into three steps: we first construct a sequence (a k ) k >i of surface tensions 
approximating a and enjoying the properties (A4)-(A6) for k > 4 (with constants depending of course on k). Owing 
to this construction, we may apply Theorem [3] to obtain, for each k > 4, a nonnegative weak solution (hk,T k ) to 
G])— (j4j) satisfying (|23|) . We then show that (h k , T k )k>4 is compact in suitable function spaces. In the last step, we 
identify the equations satisfied by the cluster points (h, T) of (hk,T k ) k >4- 

4.1. Construction of approximate surface tensions. For k > 1, we set <7fc(l) := and 

r 



ka' 



a'(s) 



a'(k) 



k l+6 



for 
for 
for 



1 

S< k' 



1 



< s < k. 



k 

s > k. 



(41) 



Recall that 9 is defined in ([8)1. Denoting a family of even mollifiers by (xe)e>o, we introduce then the approximate 
surface tension a k by 

a' k := xi/k* * a' k , <Xfc(l) = 0. (42) 
The following proposition verifies that we can apply Theorem [3] to any approximate surface tension. 

Proposition 7. Given k > 4, the free energy g k := g ak associated to o~k via formula ([7]) satisfies (A4)-(A6). 

Proof. By construction, a' k G C°°(R) and, owing to the properties of Xi/k 2 i straightforward computations yield that 



= 

400 = 



a'(k) 



for all s < 



fc-1 



k i+e 



for all s > k ■ 



1 

w 



(43) 
(44) 
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In particular, it follows from g3]) that [s ^ a' k (s)/s] G C°°(M) and g k satisfies (A4). Next, if s G ((fc-l)/fc 2 , fc + (l/fc 2 )), 
we have s - (l/k 2 ) > (k - 2)/k 2 > l/(2fc) and it follows from (0TJ), (HI), and (H2) that 



"*00 < 



i/fe,fe] 

We then infer from (|35 ]) -(|S5 |1 that 



kcr' r l( ,i/fc)(r) + cr'(r) l(i/ fe , fe )(r) + cr'(k) l( fe ,oo)(r) 

^ CT ' \%) 1 ( ^ 1 /fc)( r ) + ^ ^'W 1 (i/k,k){r) + i cr'(fe) l( fc ,oo)W 

If 1 

< - sup {<t'} / Xi/fc 2 (s - ?')dr = - sup {er'} < . 

z fi/fe,fel Jk z fi/fe,fel 



Xi/fc 2 - r) dr 

Xi/fc 2 (s- r) dr 



(45) 



k 



if s < 



jfc - 1 



> { 



1 r n -r 1 

1 -r , 1 

if k + — < s , 



k i+e k 2 

and we obtain the existence of a constant c k > for which (A5) holds. Finally, it follows from ((5]) that, for s <G 
((k-l)/k 2 ,k + (l/k 2 )), 



ffc(s) > -/ [(1 + ct ) l( ,i/fc)( r ) + CT o l(i/fc,fc)( r ) + (°"o + prpe j l(fe,oo)Wj Xi/k--(s-r)dr 



> -(2 + oib) / Xi/fc»(s-r)dr = -(2 + £r ). 

Noting that ([5} and (|4"5)) - (|4"4")) guarantee this lower bound also for s € [0, (k— l)/k 2 ) and s > k+ (l/k 2 ), we conclude 
that 

4(s) > -(2 + cr ) for s>0. (46) 
In addition, it follows from dHJ and (|43|) that a' k (s)/s > —(1 + ctq) & for s G (-co, (k — l)/k 2 ]. These two facts give 



k — 1 

fe(l + (7 ) for s < 



< 



k 2 (2 + a a ) 



for s > 



k 2 
k-1 



k-l " ~ k 2 

and we obtain (A6) with r = (so that it also holds true for arbitrary r G (0, 2)). 



(47) 



□ 



The previous proposition and Theorem [3] ensure that, for any T > and k > 4, there exists at least a nonnegative 
weak solution (hk,^k) to |T])~((4|) with surface tension a k and initial condition (ho, To). We prepare the study of 
compactness properties of the sequence (h k , Tk)k>4 by deriving technical properties of the approximate surface tensions 

(o~k)k>4- 

Proposition 8. If 9 is the exponent given by ([5)). then there exist constants G\,C<i € (0,oo) such that, for k > A, 

<g k (s) < Ci (1 + s 2 ) /or all s>0, (48) 
/v.s 



- (2 + ,o)<^)<-^ II+s)e(1 + M 
Moreover, (o~ k )k>A converges uniformly to a on compact subsets of [0, oo). 



/or a// s > . 



(49) 
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Proof. The first inequality in (|49[) having already been proved in (|46j) . we concentrate on the second inequality and 
first establish a similar estimate for a' k . As the surface tension a satisfies (f5T) , there holds 



3i(*) < 
5i(s) < 



o-i 



(l + s e ) 



< 



k i+e 



< 



(1 + s) 6 
1 



whence 



(1 + s) 6 
Cks 



for s £ 
for s > k , 



since fcs/ (1 + fcs) < 1 for s > 0. Also, 



(l + fcs)(l + s) e 
A'.s 



for s > i 
fc 



for s e 



(l + fcs)(l + s) 

Consequently, if s > (fc — l)/fc 2 , we have s — (1/fc 2 ) > (k — 2)/k 2 and, as k > 4 



0, 



we have 



, <7(fc(s-£)) C 2 fcs 

ky> ~ (1 + H± + s))(l + (± + s))e - (l + ks)(l + s) 6 

2 



for s > 



Since a' k (s) = a' k (s) for s < (k - l)/k 2 by (g3J|, we end up with 



°100 < 



C 2 fcs 



for s > , 



{l + s) e {l + ks) 

and thus obtain (|4"!)|). We next note that, given R > and s <G [0, R], it follows from © that, for k > R, we have 



|crfc(s) - er(s)| 



< 



l/A- 

mill {s,l/fc} 
1 
fc 



r — a (r) 



c/r 



2/c 



— (7 mm < s, 



I±^»+ CT (0)-<7 



2A- 



and 



\?k(s)\ < 



l+O-Q 
k 



for s £ 



Consequently, owing to the continuity of a in [0,oo) and the properties of the convolution, the sequences {&k)k and 
(cfc)fc converge uniformly to a on compact subsets of [0,oo). 

Finally, integrating the bound (|46)) gives <?fc(s) < (2 + oo) (sins — s + 1) < (2 + ao)(l + s 2 ) f° r s > 0, whence 
(H§. □ 

4.2. Compactness. Let T > 0. The main difference here with the strategy employed in Section 13.21 is that we no 
longer have an estimate on (Tk)k in £oo(0, T; L 2 (0, 1)) but only in ioo(0, T; £i(0, 1)), and this requires a different 
approach to the compactness issue for (Tk)k- Let us collect the estimates available for (hk,Tk)k which result from 
([T|)-(|4]), ((23)) . and the nonncgativity of gt- 
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(1) Conservation of matter: for t G [0,T], there holds 



h k (t, x)dx = / ho(x)dx 



(2) Energy estimate: for t € [0,T], there holds 



Tk(t, x)dx = / To(x)dx 



(50) 



T r l 



\d x h k (t,x)\ 2 dx + 



o Jo 



hl(s,x)l (0tOO )(hk{^x)) „ 2 fe fc (a,g) , n 2 
|y x ft fe (s,o;)| H |9 x o- fe (r fe (s,x))| 



21 



deeds 



T r l i 



D 



a' k (T k {s,x)) 



|a 8 r fc («,a;)| 2 dxd8< / 
Jo 



\d x h (x)f 



9k{F (x)) 



(51) 



dec. 



/o Jo r fe (s,x) 

Moreover, both /ife and Tfc are nonnegative a. e. in (0, T) x (0, 1), and ||.9/c(ro)| i < C\ ( 1 + j | To j || ) by (J3HJ) ■ Consequently, 
([ST))) and (f5"T| . together with the lower bound (|19"|) on — er^, and the Poincare inequality yield: 

(B.l) (h k ) k is bounded in L oo (0, T; H 1 ^, 1)) and (r fc ) fc is bounded in L M (0, T; Li(0, 1)). 

(B.2) {hl /2 l {0 ^ } (h k )d 3 x h k ) k , (d x r k /(l + T k )^/ 2 ) k , and {y/h~ k d x a k {T k )) k are bounded in L 2 ((0,T) x (0,1)). 
We then infer from (|20[) (with surface tension <jfc), (B.l), and the embedding of ff 1 (0, 1) in Loo(0, 1) that 

(h k )k is bounded in £^((0, T) x (0, 1)) and (d t h k ) k is bounded in L 2 {0, T; (i? 1 (0, 1))*) . (52) 

Next, we prove the following embedding: 

Lemma 9. Let T be a nonnegative function in Li(0, 1) such that (1 + r)( 1 ~ e '- ) / 2 6 H 1 ^, 1). XTien i/iere exists Co < oo 
depending only on 9 such that, after possibly redefining T on a set of measure zero, T 



||r|| c o,(i-e)/2([ 0jl ]) < Ce 



1+ [ T(x)dx 1+ / 
Jo Jo 



Proof Set 



G := 



o (i + r(.x))d+«) 
1 IWx)| 2 



ri.7- 



o (l + T(x)y+o 



dx < 



(l-#) 2 ' 

We assume T to be smooth for simplicity and focus on the distance yjl + T(x) ~ y/l + T(y) for < x < y < 1. Then, 
by Holder's inequality 



VTTW) - y/T+W) 



< 



l«z)| 



Vg 



x VT+T(zj 

n/2 



dz < 



\d x T(z) 



T dz 



1/2 



(1 + T(z)) e dz 



(l + T{z))dz 



x (i + rw+o 

\y-x\^/ 2 <VG (1 + llrllO / 2 \y-x\^ 2 . 



1/2 



Since 



v /i + r(z)dz<(i + ||r|| 1 ) 1/2 , 



integrating the above inequality with respect to y over (0, 1) ensures that ||vT+T'|j co < (1 + j|rj| i) 1 ^ 2 +v / G (1 + ||r||i) e ^ 2 , 
so that there exists Cg depending on 9 only such that 



||Vl + r|| C o,(i-e)/a([ 01 ]) < Cg 1+ T(x)dx 

Jo 

We conclude using the classical trick T = (y/l + T) 2 — 1. 



1/2 r 



\d x T(x) 



o (i + r(x)) 



1+0 



dx 



1/2 



□ 
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Now, (B.l), (B.2), and Lemma [9] yield that 

(r fc ) fc is bounded in Loo(0, T; Li(0, 1)) R Li(0, T; C°'( 1_s )/ 2 ([0, 1])) . (53) 

In particular, since 1 1 1 1 i — II^Hoo l|rfc||i, we have that 

(T k ) k is bounded in L 2 ((0,T) x (0, 1)) . (54) 

Owing to ([4"5]). a first consequence of ([EH) is that (<T k (T k )) k is also bounded in L 2 ((0,T) x (0,1)). Furthermore, it 
follows from ([ST]). (|49)l . and ([54]) that 



< 



=4P) 2/3 |a.r fc | 4 / 3 

r T dxds] I f f (T k \a' k (T k )\f dxds 

JO Jo 1 k I \Jo Jo 



[ [ \d x a k (r k )\ 4/3 dxds = [ [ ( |^r fc | 4/3 (r fc |^(r fc )|) 2/3 dxds 

Jo ^0 JO JO \ *• k J 



1/3 

< C(T) (2 + a ) 2 



J^Tldxds) <C(T) 



Consequently, 



K(r fc )) fe is bounded in L 4/3 (0, T; W^O, 1)) . (55) 

Finally, ([211) (with surface tension a k ), (B.l), (B.2), (|5^|) . and ([53")) guarantee that 

(9 t r fe ) fc is bounded in L^O, T; (#^(0, 1))*) , (56) 

where H^(Q,1) := {w £ Lf 2 (0,l) : d x w(Q) = d x w(l) = 0}. Hence, owing to the compactness of the embeddings 
of H 1 ^, 1) and e^-^QO, 1]) in C([0, 1]) and the continuity of the embedding of C([0, 1]) in either (iJ^O, 1))* or 
{H$f(0, 1))*, we infer from (B.l), (|52]l . ([55]) . ([SB"]), and [121 Corollary 4] that there are a subsequence of (h kl T k ) k (not 
relabeled) and functions h and T such that 

/!^/iinC([0,T]x[0,l]), r fc ->rinLi(0,T;C([0,l])). (57) 

In addition, (d x h k ) k being bounded in Loo(0, T; L 2 (0, 1)) by (B.l) and (d x a k (T k )) k being bounded in L 4 / 3 ((0,T) x 
(0, 1)) by (|55p . we have, up to an extraction of a subsequence and for some function S, 

d x h k ^d x h weakly-* in ^ (0,T;L 2 (0,1)) and 5 K c7 fe (r fc )^S in L 4/3 ((0, T) x (0, 1)) . (58) 

As a consequence of (B.l), ([54]l . (fSTj) . and ([58]) . we get that the limits satisfy 

h G Loo (O.Tjfl^O, 1)), fe>0, Tg Loo(0,T;L 1 (0,l))nL 2 ((0,T) x (0,1)), T>0. (59) 

Finally, thanks to (B.l), l[58]). and fS"9"]). we have 



/ |r fc (t,x)-r(t,x)| 2 dxdt < sup {||r fc (s)||i + ||r(a)||i} 

o Jo se[o,T] 



A;— >oo 



so that we also have 

r fe ^rinL 2 ((0,T)x (0,1)). (60) 



THIN FILM MODEL WITH CAPILLARY EFFECTS AND INSOLUBLE SURFACTANT 



15 



4.3. Identifying the limit system. According to the uniform bounds (B.l), (B.2), and ([52]) . we first obtain that, 
up to an extraction of a subsequence and for some function j 1 and j 2 , 

ftfc /2 l( ,oo)(/i*)^fc^Ji, hl /2 d x a k (T k )^: h , inL 2 ((0,T)x (0,1)). (61) 

Arguing as in [3J Section 3] and [7J Section 3.4], we first deduce from (B.2) and ([57]) that d^h belongs to L 2 ('P(5)) for 
all S > where 7>(<J) := {(*, a;) € (0, T) x (0, 1) : h(t, x) > 5} and j x = h s / 2 d*h in {(i, i) G (0, T) x (0, 1) : x) > 0}. 
Combining this result with f|57"[) yields 

Urn / / h 3 k l (0 ,oo)(hk)d^h k d x (dxds = [ [ h 3 l {0iOo) (h) d%h d x £dxds , 

k^Jo Jo Jo JO 

lim / / hll { o !0o) {h k )T h dlh h d x (dxds= I { h 2 l {Qt00) {h)Td 3 x h d x (dxds, 

for any ( G C°°([0, T] x [0, 1]) such that ((T, x) = for all a; G [0, 1] and 9 X C(*. x) = for all (t, x) G [0, T] x {0, 1}. We 
next claim the strong convergence 

<7fc(r fc ) ^ tf(r) inL 2 ((0,T)x (0,1)). (62) 
Indeed, on the one hand, we readily infer from (|49)l and ((60)) that 



/ / K(r fc )-a fc (r)| 2 dxdt < (2 + a ) 2 [ [ |r fe -r| 2 da;di —t 0. 
Jo Jo Jo Jo fe ^°° 

On the other hand, it follows from Proposition [S] and (g5J) that <T fe (r) -> <r(r) a.e. in (0,T) x (0,1) with |cr fc (r)| < 
(2 + (To) (1 + r) G L 2 ((0, T) x (0, 1)), whence 

lim / / |cr fc (r) - <r(r)| 2 dxdt = 
fc ^°°Jo Jo 



by the Lebesgue dominated convergence theorem. Thus, (|62|) holds true. 

Together with ([55]). the convergence (JB2J) ensures that cr(r) G i 4/3 (0, T; ^^(0, 1)) and E = 9 x cr(r) in 

Next, collecting d57j, (JSHJ) , and (pj yields j 2 = ft^^cr^), so that ft 1 / 2 ^^^) G L 2 ((0,T) x (0,1)). It is then 
straightforward to pass to the limit as k — > oo in the remaining terms in the weak formulation (|20 |) -(|21 j) for (hk,Tk) 
and conclude that (h,T) is a weak solution to ([I])~(j4]) with surface tension a and initial data (ho, To). This completes 
the proof of Theorem [2] 

Remark 10. We shall point out that our strategy to prove Theorem [2] by approximating the surface tension a G 
C([0, oo)) (~1 C 1 (0, oo) 6?/ surface tensions (o-k)kGN satisfying the assumptions of Theorem [3] does no£ yield existence 
of nonnegative weak solutions for the limiting case 6 = 1. More precisely, if 8 = 1, then the analogue of Lemma [9] 
for the energy dissipation merely yields a control on %/T in the space of continuous functions (instead of a Holder 
space as in the case 8 G [0, 1) ), and we thus lose compactness of the concentration of any bounded family of solutions 
in L 2 ((0,T) x (0,1)). It seems that this threshold is of high importance. Indeed, provided — o~'(T) is dominated by 
1/(1 + r) at infinity, a good choice of multiplier for ([l])-([2]) yields that the integral (|24j) (with 6 = 1) measures the 
dissipation of energy for any small solution. 
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